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Abstract. Low-mass X-ray binaries, recycled pulsars, cataclysmic variables and magnetically 
active binaries are observed as X-ray sources in globular clusters. We discuss the classification 
of these systems, and find that some presumed active binaries are brighter than expected. We 
discuss a new statistical method to determine from observations how the formation of X-ray 
sources depends on the number of stellar encounters and/or on the cluster mass. We show 
that cluster mass is not a proxy for the encounter number, and that optical identifications are 
essential in proving the presence of primordial binaries among the low- luminosity X-ray sources. 
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1. Introduction 

The first celestial maps in X-rays, in the early 1970s, show that globular clusters 
harbour more X-ray sources than one would expect from their mass. As a solution to 
this puzzle it was suggested that these bright {Lx lO'^^erg/s) X-ray sources, binaries 
in which a neutron star captures mass from a companion star, are formed in close stellar 
encounters. A neutron star can be caught by a companion in a tidal capture, or it can 
take the place of a star in a pre-existing binary in an exchange encounter. Verbunt & Hut 
(1987) showed that the probability of a cluster to harbour a bright X-ray source indeed 
scales with the number of stellar encounters occurring in it; whereas a scaling with mass 
does not explain the observations. 

With the Einstein satellite a dozen less luminous {L^ ^ 10'^^ erg/s) X-ray sources were 
discovered in the early 1980. ROSAT enlarged this number to some 55, and now thanks 
to Chandra we know hundreds of dim X-ray sources in globular clusters. The nature and 
origin of these dim sources is varied. Those containing neutron stars, i.e. the quiescent 
low-mass X-ray binaries in which a neutron star accretes mass from its companion at 
a low rate and the recycled or millisecond radio pulsars, have all formed in processes 
involving close stellar encounters. The magnetically active binaries, on the other hand, 
are most likely primordial binaries, with stars that are kept in rapid rotation via tidal 
interaction. Cataclysmic variables are binaries in which a white dwarf accretes matter 
from a companion. In globular clusters they may arise either via stellar encounters, or 
from primordial binaries through ordinary binary evolution - this is expected to depend 
on the mass and density of the globular cluster. 

In this paper we describe the classification and identification of the dim sources in Sec- 
tion 2, and make some remarks on the theory of their formation in Section 3. In Section 4 
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Figure 1. Left: X-ray hardness-luminosity diagram for dim sources in globular clusters. I: qui- 
escent low-mass X-ray binaries, II: cataclysmic variables III: cataclysmic variables and magneti- 
cally active binaries. From Pooley & Hut (2006). Right: Colour-magnitude diagram of NGC 6752 
on the basis of HST-WFPC2 data; objects within X-ray position error circles are marked. Left of 
the main sequence we find cataclysmic variables, above it active binaries. Updated from Pooley 
et al. (2002a). 

we will discuss a new, and in our view more accurate, way to compare the numbers of 
these sources with theoretical predictions. 



2. Clasification and identification 

Work on the dim sources is progressing along various lines. Grindlay and coworkers 
study one cluster, 47Tuc, in great detail (Grindlay et al. 2001, Edmonds et al. 2003, 
Heinke et al. 2005). Webb and coworkers use XMM to obtain high-quality X-ray spectra 
(e.g. Webb et al. 2006, Servillat this meeting). Dim sources are also found in clusters in 
which individual sources are the main target, such as Terzan 1 and 5, and M 28 (Wijnands 
et al. 2002, Heinke et al. 2003, Becker et al. 2003). Lewin initiated a large program to 
observe clusters with very different central densities and core radii, and thereby to provide 
material for tests on the dependence on these properties of the numbers of dim sources. 
Further references to all this work may be found in the review by Verbunt & Lewin 
(2006); and in the remainder of this Section. 

The first classification of the dim sources may be made on the basis of the X-ray prop- 
erties only (Fig.[T]). The brightest sources in the 0.5-2.5 keV band, at 10"^^ erg/s, 
tend to be quiescent low-mass X-ray binaries. To better use the Chandra range, one 
may also select the brightest sources in the 0.5-6.0 keV band, and select soft sources, 
with a high ratio of fluxes below and above e.g. 2keV: /o.5-2.okov//2.o-6.okcV 1- 
Such sources also are mostly quiescent low-mass X-ray binaries. Between 10'^"'^ and 10'^^ 
erg/s most sources are cataclysmic variables, especially when they have hard spectra 
/o.5-2.okcv//2.o-6.okcV < 1- Thc faintest sources include magnetically active binaries, of- 
ten with soft X-ray spectra. For many faint sources the number of counts is too low to 
decide on the hardness of the spectrum. 

The second step in classification can be made when identification with a source at 
otiier wavelengtiis is made. Positional coincidence of an X-ray source with the accurate 
radio position of a millisecond pulsar provides a reliable identification and classification. 
Positional coincidence with optical sources is only significant if the highest possible as- 
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Figure 2. Left: X-ray luminosity as a function of absolute visual magnitude for nearby stars 
{selected from Hiinsch et al. 1999, for details see Verbunt 2001) and for RS CVn systems (from 
Dempsey et al. 1993). The upper bound Eg. 12. 31 is indicated with a solid line. We convert the 
X-ray fluxes in the 0.1-2.4 keV range (scale on the left) to the 0.5-2.5 keV range by multiplication 
with 0.4 (scale on the right). Right: X-ray luminosity as a function of absolute visual magnitude, 
for dim X-ray sources in globular clusters. The assumed separatrices Eqs. 12.112.21 are indicated 
with dotted lines, the upper bound Eg. 12. 31 with a solid line. It is seen that some X-ray sources 
classified as active binaries in globular clusters are well above this bound. 

trometric accuracy is used to limit the number of possible counterparts (e.g. Bassa et 
al. 2004). The position of these possible counterparts in a colour-magnitude diagram is 
then used to select the probable counterparts. Cataclysmic variables are bluer than the 
main sequence stars, and magnetically active binaries may lie above the main sequence. 
Systems on the main sequence cannot be unambiguously classified: they may either be 
cataclysmic variables in which the optical flux is dominated by the donor star, or main- 
sequence binaries with unequal masses whose optical light is dominated by the brighter 
star. If a periodicity is found in the X-rays that corresponds to a period at another wave- 
length, e.g. the pulse period of a pulsar or the orbital period of a binary, identification 
and clasification are secured simultaneously (e.g. Ferraro et al. 2001). 

A very useful discriminant in X-ray astronomy in general is the X-ray to optical flux 
ratio. In the case of globular clusters we can use the known distance to determine the 
optical to X-ray luminosity ratio (Fig[2]). On the basis of in particular the extensive 
data on 47Tuc (Edmonds et al. 2003 and references therein) one finds that the lines of 
constant optical to X-ray luminosity ratio 

log Lo.5-2.5kcv( erg/s) = 36.2 - 0.4Afv (2.1) 

separates the quiescent low-mass X-ray binaries above it from the cataclysmic variables 
below. The line 

logLo.5^2.5kcv(erg/s) = 34.0 - 0.4Mv (2.2) 

roughly separates the cataclysmic variables from the magnetically active binaries. 

This latter separatrix leads to a surprise when one compares it with the X-ray lumi- 
nosities of nearby stars and of known magnetically active binaries, i.e. RS CVn systems, 
near the Sun. For main-sequence stars in the solar neighbourhood, the X-ray luminosity 
increases with the rotation speed, up to an upper bound given approximately by 

log Lo.5-2.5kcv( erg/s) - 32.3 - 0.27Mv (2.3) 

as illustrated in Fig. [5] (left). This bound is lower than the separatrix given by Eg. 12. 21 
especially for brighter stars. This would imply that active binaries in globular clusters 
can have higher X-ray luminosities than similar binaries near the Sun. We suggest, how- 
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ever, that the classification must be reinvestigated, and that some of these objects are 
cataclysmic variables. The absence of the blue colour expected for a cataclysmic variable 
(see Fig.[T|) then requires explanation - e.g. as a consequence of the non-simultaneous 
measurements at different colours combined with source vari ability^ 



3. Some remarks on theory 

Binaries in a globular clusters change due to their internal evolution and/or due to 
external encounters. To describe the current cluster binary population one must track 
the events for each primordial binary and for each binary that is newly formed via tidal 
capture, throughout the cluster. The first estimates of the formation of binaries with 
a neutron star necessarily made a number of drastic simplifications. The sum of all 
encounters (of a neutron star with a single star, or with a binary) was replaced with an 
integral over the cluster volume of the encounter rate per unit volume. Four assumptions 
followed: the number density ni, n2 of each participant in the encounter scales with the 
total mass density p, the relative velocity between the encounter participants scales with 
the velocity dispersion w, the interaction cross section A is dominated by gravitational 
focussing so that A cx and the encounter rate is dominated by the encounters in 

the dense cluster core. Hence one writes the cluster encounter rate F' as 

F' = / nin2AvdV cx Po^r^ jv (3.1) 
Jv 

where Po is the central density and the core radius. If one further eliminates the 
velocity dispersion through the virial theorem, v cx Tc-y^p^, one has 

F' cx po^-^c' = r (3.2) 

where F is referred to as the collision number. With a life time r the expected number 
of binaries of a given type is 

TV = F'r cx Ft (3.3) 

A major advantage of these simple estimates is the clear connection between (the 
uncertainty in) the input and (the uncertainty of) the output. Thus, if 7ii,rt2 are the 
number densities of neutron stars and of binaries, respectively, Eqs. l3.1ll3.3l indicate that 
the uncertainty in the number N of neutron star binaries scales directly with the uncer- 
tainties in n\ and n2. Similarly, if we overestimate the life time r of a binary by a factor 
10, the estimated number N is overestimated by the same factor. 

Thanks to a concerted effort by various groups fairly detailed computations of the 
happenings in globular clusters are now undertaken. This is a fortunate and necessary 
development, as many details cannot be understood from the simple scalings above. For 
example, the wide progenitors of cataclysmic variables are destroyed by close encounters 
before they evolve in a dense cluster core (Davies 1997), but evolve undisturbed into 
cataclysmic variables in the outer cluster regions, from where they can sink to the dense 
core to form a significant part of the current population there (Ivanova et al. 2006). 

A disadvantage of complex computations is that they tend to hide the uncertainties. 
If the cross sections A are described in paper I of a series, and the life times r in paper 
III, the large uncertainties in them tend to be less than obvious in paper V where the 

f In his contribution to this meeting. Christian Knigge shows that the tentative counterpart 
of W24 in 47 Tug, a possible active binary according to Edmonds et al. (2003, Sect. 4. 5), has blue 
FUV-U colours, which suggests that it is a cataclysmic variable. At My = 2.6 and = 8.7 x 10^° 
erg/s, it actually lies below the line given by Eg. 12. 31 
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Figure 3. Optical (*, scale on left) and 
X-ray (+, scale on right) lightcurves of the 
dwarf nova YZ Cnc through several outburst 
cycles: the outburst in the optical luminos- 
ity is accompanied by a marked drop in the 
X-ray luminosity. After Verbunt et al. (1999) 

final computations are described. The confidence expressed in summaries of the results 
of such computations is sometimes rather larger than warranted. An uncertainty in A 
or r has an equally large effect in complex computations as in simple estimates. As a 
further illustration we discuss two other uncertainties. 

The first relates to the question what happens when mass transfer from a giant to its 
binary companion is dynamically unstable. It is usually assumed that a spiral-in follows, 
in which the companion enters the envelope of the giant and expells it through friction. 
The outcome of this process is computed using conservation of energy, which implies a 
drastic shrinking of the orbit (Webbink 1984). However, the study of nearby binaries 
consisting of two white dwarfs shows that the mass ratios in them are close to unity 
(e.g. Maxted et al. 2002). Such binaries can only be explained if the consequences of 
dynamically unstable mass transfer are governed by conservation of angular momentum, 
rather than by the energy equation. If the mass leaving the binary has roughly the 
same specific angular momentum as the binary, the orbital period changes relatively 
little during the unstable mass transfer and concomitant mass loss from the binary (Van 
der Sluys et al. 2006). The standard prescription of dynamically unstable mass transfer 
hitherto implemented in globular cluster computations must be replaced. 

The second uncertainty relates to the conversion of the mass transfer rate in a cata- 
clysmic variable to the X-ray luminosity. This conversion does not affect the evolution 
of the binary but it is important for comparison with observations, as most cataclysmic 
variables in globular clusters are discovered as X-ray sources. It is generally assumed that 
the X-ray luminosity scales directly with the mass transfer rate: Lx oc M. Alas, reality 
is more complicated, and indeed in most cases the X-ray luminosity goes down when the 
mass transfer rate goes up. This is demonstrated unequivocally in dwarf novae whose X- 
ray luminosity drops precipitously during outbursts (Fig.|3]), but there is evidence that it 
is true in the more stable nova- like variables as well (Verbunt et al. 1997). But exceptions 
are also known: the dwarf nova SS Cyg has higher X-ray flux during outburst than in 
quiescence (Ponman et al. 1995). Theoretical predictions of the numbers of cataclysmic 
variables in globular clusters that radiate detectable X-ray fluxes, are not believable when 
based on proportionality of X-ray flux and mass-transfer rate. 

Observational evidence for the numbers of binaries of various types, and of the depen- 
dence of these numbers on cluster properties, may be collected and used to constrain the 
theories on formation and evolution of various types of binaries in globular clusters (e.g. 
Pooley et al. 2003, Heinke et al. 2006, Pooley & Hut 2006). 

In the next Section we describe a new, and we hope more accurate, method of analysing 
source numbers: this method is based on direct application of Poisson statistics. This topic 
brings one of us, FV, to a brief Intermezzo. 
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MEGfiLiTHiC 
SCIENCE 



DOUGLAS C. HEGGIE 



Ancient Mathematics and Astronomy 
in Northwest Europe 




Figure 4. Left: the classic book by Douglas Heggie on Megalithic Science. Right: Sketch of 
Stonehenge with the circular bank & ditch (the henge) , the rectangle of the four Station Stones 
within it, and the directions towards the extreme risings of Sun and Moon at the winter and 
summer solstitia of 2500 B.C. and 2000 A.D. The minute shift in 4500 yr is due to the small 
change in obliquity. 

Intermezzo by FV: another side of Douglas Heggie 

The first time that I had need of understanding of binomial and Poisson statistics was in 
the early 1980s, when I bought a book by Douglas Heggie (1981) called Megalithic Science, 
Ancient Mathemathics and Astronomy in Northwest Europe. The first half of this book 
discusses the question whether the megalith builders used a standard measure of length 
(the answer is no), and the second half studies the question whether megalithic structures 
had astronomical orientations. Whereas Douglas is still rather modest about this book, it 
is still the best introduction to the issues in megalithic astronomy (McCluskey 1998, p. 11). 
It, for the first time, explained the statistics of the study of astronomical orientations 
of megalithic structures. In my view it is not an exaggeration to state that archaeo- 
astronomy may be divided in a pre-Heggie era, in which the statistics is usually wrong, 
and a post-Heggie era, in which at least some people, following Douglas's prescription, 
do their statistics right. Among the latter I may specifically mention Ruggles (1999). 

As an example of an application of binomial statistics I paraphrase the discussion on 
p. 198 of the book, of the four Station Stones in Stonehenge. These Stones were placed, 
near 2500 B.C., in a rectangle with the short sides parallel to the main axis of Stonehenge 
(Fig. [3]). This main axis points from the center of the rectangle through a gap in the 
bank-plus-ditch (the henge) surrounding it toward the location on the horizon of the 
most northern annual rising of the Sun, at summer solstice. The Moon ranges in a 19 yr 
cycle from about 5° north of the Sun to 5° south of the Sun, which leads to a range of its 
most northern and southern annual risings. The extremes of these ranges, together with 
the two solstitia, define 6 positions on the eastern horizon. Accepting a range of ±1.8° 
for each position, we find a probability that an arbitrary direction towards the eastern 
horizon hits one of these positions p = 6 x 2 x 1.8/180 = 0.12. The rectangle of the Station 
Stones defines three new directions towards the East: along the long side and along the 
diagnonals. Two of these hit one of the 6 positions: extremes of the southernmost annual 
setting of the Moon in the 19 yr cycle. The probability of 1 chance hit in 3 trials with 
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p — 0.12 is 28%, the probability of 2 chance hits 4%. Most hkely one hit is due to chance, 
and one hit intentionaL 

Thanks to new research at Stonchenge, done after the pubhcation of Megalithic Sci- 
ence, we now know that the central structure, circle and horse-shoe form, were put up 
simultaneous with the Station Stones, which implies that the view along the diagonals 
was blocked. The intentional hit is therefore the one along the long side. Since the long 
side of a rectangle is necessarily perpendicular to the short side, this implies that the 
location of Stonehenge was selected for its latitude. 

4. Testing models against observations 

For the denser clusters it is found that the X-ray sources lie well within the half- 
mass radius (e.g. NGC 6440, Fig.l of Pooley et al. 2002b), whereas for (apparently) large 
clusters only the area within the half-mass radius is covered (e.g. 47 Tuc, Grindlay et al. 
2001). Thus, the analysis generally deals with the sources within the half-mass radius. 
If Nfi is the observed number of X-ray sources, and Nf, the number of background (or 
foreground) sources not related to the cluster, the number of cluster sources within the 
half-mass radius is Nc = Nh — Ni,. From a model we may obtain an estimate of the 
expected number of cluster sources /ic, and we also may estimate the expected number 
of background sources fit- 

In clusters with a small number of sources, Nc ^ 10, say, one cannot apply chi-squared 
statistics. One way of solving this is by adding such clusters together, to obtain sufficiently 
large numbers. This is done by Pooley & Hut (2006). In this process, the information on 
the indicidual clusters is lost. 

To avoid this problem, we fit as follows. The probability of observing N when /i is 
expected according to a Poisson distribution is 

P{N,f,)^^e->^ (4.1) 

An important aspect of the Poison function for our application is its asymmetry: P(3, 0.1) ^ 
P(0, 3), i.e. the probability of observing 3 sources when 0.1 is expected is very much 
smaller than the probability of observing sources when 3 are expected. We now con- 
sider 

Nh = Nc + Nb = ^ic + IJ-b with Hc = ar + bM (4.2) 
Both Nc,Ni, are realizations of Poisson distributions given by /Xc and fib respectively. Hb 
is determined for each cluster separately, e.g. from the observed number of sources well 
outside the half-mass radius. For each value of Nc we take Nb = Nh — Nc, and then 
compute the combined probability P{Nc, fic)P{Nb, fib)- We then select the Nc,Nb pair 
with the highest combined probability (Fig. [5]). The fitting procedure consists of varying 
a and b to maximize 

P = l[[P{Nc,fic)P{Nb,fib)]j (4.3) 

i 

where j indexes the clusters. 

Table[l] lists the numbers that we use in our fitting. We note from the Table that mass 
cannot be used as a proxy for collision number. The cluster with the highest collision 
number has a ratio F/M ~ 7 (normalized on M4), the cluster with the highest mass has 
T/M ~ 0.2. This suggests that mass and collision number can be well discriminated, and 
we confirm this below. The derived values for Nh and fib (for sources with io.5_6.okcV > 
4 X 10^° erg/s) correlate with the derived values of F and M, since all are based on 




Figure 5. Three examples of the combined fit of the number of cluster sources Nc and back- 
ground sources Ni, for one cluster. The expected number of cluster sources according to the 
model Nc = l.Sr, and the expected number of background sources fit arc indicated with each 
frame. For each realization Nc ^ Nh the accompanying realization is Nb = Nh — Nc- E.g. in 
NGC6266 the observed number Nh = 51, allowing the combined realizations of Nc,Nb as 51,0 
or 50,1 or 49,2 etc. The probability P{Nc,Hc) is indicated with a dashed line; the probability 
P{Nb,Hb) with a dotted line; and the combined probability - the product of these for each 
allowed pair - with a solid line. The cases shown are from left to right for low, roughly equal, 
and dominant background {fib .'V^, Hb ~ 0.5Nh, Hb 3> yUc). 
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Table 1. Collision number F and mass M, normal- 
ized on the values for NGC 6121, for clusters stud- 
ied with Chandra, together with the number Nh of 
sources within the half-mass radius, the expected 
contribution of background sources fib, and the num- 
ber of secure members Ns. Most values Nh and fib 
are from Pooley et al. (2003); those for NGC 288, 
NGC 6366 and NGC 6809 are from estimated conver- 
sion of 1/0.5-2.5 kcv as given in Kong et al. (2006) and 
Bassa ct al. (2007) to I/o.5-6.0kcV- Collision numbers 
and masses are derived from the central density (in 
1/0 pc^"*), core radii, and absolute magnitudes given 
in the Feb 2003 version of the Harris (1996) compi- 
lation, and they are scaled on the values for M4 = 
NGC 6121. Thus we assume that the mass to light 
ratio is the same for all clusters. 




"log P 



Figure 6. Left: distribution of " log P for 1000 random realizations of the best model fic = 1.5F, 
with the 1- and 2-a ranges indicates by solid and dashed lines. The " log P value for the best fit 
is indicated with a dotted line. Right: The best model for fic = l.ST is indicated with a solid 
line, the best number Nc with o (< indicates Nc = 0), and the observed number Nh = Nc + Nb 
with •. 
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Figure 7. 1,2,3-a contours for the best values 
for the model of Eg . 14. 21 Dashed contours are 
for the fits which allow all observed sources to 
be background sources, solid contours for the 
0.5 1 1.5 2 fits where a minimum of secure cluster sources 

a is imposed. 

assumed values for the cluster distance d and interstellar absorption Ay . For the clusters 
in the Table the correlations thus introduced are small compared to the ranges of F and 
M, because the cluster distances are relatively well known, and we ignore them in what 
follows. 

The models we fit all have the form given by Eq. l4.2l We first fit models with 6 = 0, i.e. 
assuming that the number of cluster sources depends only on the collision number. The 
best solution has a = 1.5 and is shown in Figure[6l The total probability as defined in 
Eg. 14. 31 is P — e~^^. To see whether this is acceptable, we have used a random generator 
to produce 1000 realizations by drawing random numbers from Poisson distributions with 
values /ic as given by the best model and fit as listed in Table[Tl and for each realization 
computed the total probability P. The cumulative distribution of "^logP is also shown 
in Figured the best solution is within the 95% range around the median value, and thus 
acceptable. 

Our next fit allows non-zero values for a and 6, i.e. allows a linear dependence both 
on mass and on collision number. The best solution now has a = 1.3 and b = 0.7; 
contours of 1,2,3 a around these best values are shown in FigureEl where we assume that 
the distribution of — 2('^logP — logPmax) is given by a statistic with one degree of 
freedom. (Pmax is the total probability of the best solution.) We see that the solution with 
6 = is marginally acceptable, i.e. the evidence for the mass dependence is marginal. 
Because we find an acceptable solution for F as given by Eg. 13. 2) there is no need for a 
different dependence on central density. 

Our final fit uses more information, viz. the number of (almost) certain cluster members 
among the X-ray sources, as determined though optical identifications. The importance 
of this may be seen from Fig.[5]for NGC 288. With an estimated background /.tb = 8 and, 
for a model depending on collision number only, /ic ~ 0.1, the most probable solution is 
that all 11 sources within the half-mass radius are background sources (cf. Fig. [6]). If we 
add the constraint that at least two sources included in are cluster members (Kong 
et al. 2006) this solution is no longer possible, and the post probable solution now has 
Nc — 2 and Ni, — 9. Because the collision number is so small, any cluster source must 
be a primordial binary, as indeed argued by Kong et al. The best solution for all clusters 
combined now has a = 1.2 and b = 1.1, and the solution with 6 = now lies well outside 
the 3-(T contour, i.e. the dependence on mass M is significant (Fig.[7])- 

5. Conclusions 

• mass M is not a proxy for collision number F 

• the number of dim sources scales both with collision number F and with mass M 
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• scaling with mass only is not acceptable 

• correct treatment of the background is important, esp. for faint sources 

• to prove the mass-dependence optical identifications are essential 
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